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Abstract 



The modified associate formalism is compared with similar models, such as the classical associate model, the associate species model and the 
modified quasichemical model. Advantages of the modified associate formalism are demonstrated. 
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1. Introduction 
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In this paper we compare the modified associate formalism 
(MAF) proposed in the first paper [1] with the classical asso- 
ciate model (CAM) described by Prigogine and Defay [2] and 
the associate species model (ASM) described by Besmann and 
Spear [3]. Where it is possible, we compare MAF with the 
classical quasichemical model (CQM) of Fowler and Guggen- 
heim [4]. We also compare MAF with two modifications of 
CQM. The first modification suggested by Pelton and Blander 
[5] include: 

selection of the coordination numbers that allow one to set 
the composition of maximum ordering so as to comply with 
the experimental data; 

- an empirical expansion of the molar Gibbs energy of the qua- 
sichemical reaction as polynomials in terms of coordination 
equivalent mole fractions of solution components, which is 
used for fitting available experimental data. 

The first version of the modified quasichemical model is re- 
ferred to as MQM(1986) in the present study. Further modifi- 
cations of the quasichemical model introduced by Pelton et al. 
[6] are: 

- expansion of the Gibbs energy of the quasichemical reac- 
tion as polynomials in pair fractions (instead of equivalent 
mole fractions of solution components), which provides ad- 
vantages in data-fitting; 
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- dependence of coordination numbers of the solution compo- 
nents on the mole numbers of different pairs that allow one 
to set the composition of maximum ordering so as to comply 
with the experimental data for each binary system individu- 
ally. 

We refer to the second version of the modified quasichemical 
model as MQM(2000). 



2. Configurational Entropy of Mixing in Binary Solutions 

In this section, we compare the molar configurational en- 
tropies of mixing given by the solution models mentioned above 
with theoretical boundaries for the configurational entropy of 
mixing. These boundaries are discussed immediately below. 

The configurational entropy of mixing is a measure of disor- 
der and has its maximum value for an ideal solution which is 
completely disordered. For any solid or liquid solution A — B, 
the configurational entropy of mixing can not exceed that of 
ideal solution. On the other hand, when vacancies are not taken 
into account (which is the case in all considered models), the 
configurational entropy of mixing can not be negative. The val- 
ues of molar configurational entropy of mixing which are out- 
side of these boundaries are paradoxical. 

For each model, we consider the following three limiting 
cases. 

(i) The case of ideal solution of A-particles and B-particles. 
In this case, the Gibbs free energy changes on forming 
associates that consist of both A-particles and B-particles. 
The configurational entropy of mixing should reduce to 
that of the ideal solution model. 

(ii) The case of immiscible components A and B. In this case, 



the configurational entropy of mixing should equals to 
zero. 

(iii) The case of the highly ordered solution A — B. In this 
case, the Gibbs free energy change on forming of the as- 
sociates of a particular composition and of a particular 
spacial arrangement of particles approaches — °°. If the 
quasichemecal model is considered, the Gibbs free en- 
ergy change on forming (A — Z?)-bonds approaches — °o. 
The configurational entropy of mixing should be equal 
to zero at the composition of maximum ordering. 

2.1. Classical Associate Model 

Consider a solution of xa moles of the component A and 
xg = 1 — xa moles of the component B. The molar configura- 
tional entropy of mixing s con f given by the classical associate 
model [2] is expressed as 

*corf=--R(raAiM*Ai)+««iM*B 1 )+ Y, n AiBjM x AiBj))- 

(1) 

Here, Ha v nB x and riAfij are the mole numbers (in one mole 
of the solution A — B), while xa [7 *b, and xa,b ; are the mo- 
lar fractions of Ai-monoparticles, Z?i-monoparticles and A,\B ; - 
associates, respectively; R is the universal gas constant. Note 
that compositions of associates and sizes of associates (i + j) 
are determined by a modeller. The equilibrium values of the 
mole numbers are determined by minimisation of the Gibbs 
free energy of the solution subject to the mass balance and 
nonnegativity constraints (see, for example, the monograph by 
Prigogine and Defay [2] for more details). 

Now consider the limit when the Gibbs free energy changes 
on forming the associates equal to zero. As pointed out by 
Luck et al. [7], the classical associate model gives paradoxical 
values of the molar configurational entropy of mixing, which 
are larger than those given by the ideal solution model. Fig. 1, 
for example, shows the configurational entropy curves given 
by the classical associate model where only A2_B-associates 
are considered (solid line) and where A2_B-associates and AB2- 
associates are considered (dashed line). 

As shown in the figure, the curves are outside of the the- 
oretically possible range. Note also that the entropy curve is 
not symmetric with respect to vertical line with Xb= 1/2 when 
only A2_B-associates are considered. To obtain a symmetric en- 
tropy curve (which is expected in the limit of ideal solution) 
one should consider both A2_B-associates and AZ?2-associates. 

Consider the limit when the Gibbs free energy changes on 
forming all associates approach +°°. In this case, no associates 
form and the configurational entropy of mixing equals to that of 
the ideal solution model. In theory, however, the configurational 
entropy of mixing should be zero for any xg- 

Now consider the limit of highly ordered solution and assume 
that the Gibbs free energy change of forming A2B-associates 
approaches — 00. The entropy curve given by the classical as- 
sociate model in this case is shown in Fig. 2. As shown in this 
figure, the entire entropy curve is within the theoretically possi- 
ble range. Note also that the configurational entropy of mixing 
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Fig. 1 . Molar configurational entropy of mixing given by the classical associate 
model in the limit of ideal solution of A-particles and B-particles. 
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Fig. 2. Molar configurational entropy of mixing given by the classical associate 
model in the limit of the highly ordered solution. 

equals to zero at the composition of maximum ordering, since 
different spatial arrangements of particles in an A2_B-associate 
are not considered in the classical associate model. 



2.2. Associate Species Model 

According to the associate species model [3], the configura- 
tional entropy of mixing of xa moles of the component A and 
xg = 1 — xa moles of the component B is given by 
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Here, q is the size of species (the number of particles per an 
associate specie); iiA a . q B bq are the mole numbers and XA a . q B b . q 
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are the molar fractions of the A^B^-species. Similar to the 
classical associate model, the size and compositions of species 
are determined by a modeller. 

Consider the limit of ideal solution, assuming that the so- 
lution consists of the associate species A q , B q , A 2q pB q p and 
A q pB 2q p. The Fig. 3 shows the entropy curves given by the 
associate species model for q= 1, q = 2 and q = 3. In contrast 
to the classical associate model, one can adjust the configura- 
tional entropy of mixing, using the additional parameter q. For 
this particular set of species, the selection of q = 2 result in the 
entropy curve which is closed (though located outside the theo- 
retically possible range) to the entropy curve given by the ideal 
solution model. As shown in Fig. 3, the selection q = 1 results 
in substantial overestimation of the configurational entropy of 
mixing, while the selection q = 3 gives underestimated values. 

Consider the limit of immiscible components, that is, the 
Gibbs free energy changes on forming the species A 2q pB q p 
and A q pB 2q /i approach +°°. The entropy curves given by the 
associate species model in this case are shown in Fig. 4. As 
shown in this figure, the associate species model overestimate 
the configurational entropy of mixing in the limit of immiscible 
components. However, the entropy curve tends to theoretically 
expected one (s C onf = 0), when q increases. 

Now consider the limit of highly ordered solution. Assume 
also that the Gibbs free energy change on forming the associate 
species A 2q /^B q p approaches — °°. Fig. 5 shows the entropy 
curve given by the associate species model in this case. As seen 
from Fig. 5, the entropy curve with q = 3 is entirely located in 
the theoretically possible range. If q — 2, however, small part 
of the entropy curve (approximately, for < x% < 0.075) is 
located outside the theoretically possible range. The selection 
q = 1 results in paradoxical values of configurational entropy 
approximately for < xb < 0.25 and 0.6 <xb < T 




Fig. 5. Molar configurational entropy of mixing given by the associate species 
model in the limit of highly ordered solution. 

2.3. Modified Associate Formalism 

In the modified associate formalism the configurational en- 
tropy of mixing is given by (see Eq. (44) in Ref. [1]) 
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where all the notations have the same meaning as in Ref. [1], 
In contrast to previous modifications of the associate model, 
the compositions of associates that have to be considered are 
defined by the size of associates. 

Consider the limit of ideal solution. In this limit, the Gibbs 
free energies of the reactions of forming different associates 
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Fig. 6. Molar configurational entropy of mixing given by the modified asso- 
ciate formalism in the limit of immiscible components. 



Fig. 7. Molar configurational entropy of mixing given by the modified asso- 
ciate formalism in the limit of highly ordered solution. 



equal to zero. As demonstrated in the previous paper [1], the 
configurational entropy of mixing correctly reduces to that of 
the ideal solution model for arbitrary size of associates. 

Now consider the limit of immiscible components, where 
the Gibbs free energies of the associate that consist of both 
A-particles and B-particles approach +°°. In this case, Eq. (3) 
reduces to 



s conf = (xa In (x A ) + n B In (x B ) ) ■ 

m 



(4) 



where m is the size of associates. The entropy curves for m = 3 
and m = 6 are shown in Fig. 6. As seen from this figure, the 
entropy curves are entirely located in the theoretically possible 
range and tends to the theoretically expected curve (s C onf = 0) 
with increasing the size of associates. 

Considering the limit of highly ordered solution, we assume 
also that the Gibbs free energy of associate with the composi- 
tion A^B and with a particular spatial arrangement of particles 
approaches — 00. The entropy curves for m = 3 and m = 6 given 
by the modified associate formalism in this case is presented in 
Fig. 7. In this case, the curves are also entirely located in the 
theoretically possible range. 



2.4. Modified Quasichemical Model (1986) 

In the modified quasichemical model, the molar configura- 
tional entropy of mixing is given by (see Eq. (10) in Ref. [6]) 



•Sconf = -R(x A \n(x A ) + x B ln(x B )) 
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Here, xa and x B are the molar fractions of the components A 
and B; tin and xij = A,B) are the mole number and the 
molar fraction of (z — /')-pairs, respecively. The coordination- 



equivalent fractions and y B are defined as (see Eq. (6) in 
Ref. [6]) 



y A =Z A nA/(Z A nA+Z B n B ) 
y B =Z B n B /(Z A n A +Z B n B ) , 



(6) 



where Za and Z B are the coordination numbers of the solution 
components A and B, respectively. As demonstrated by Pel- 
ton and Blander [5] for constant coordination numbers and by 
Pelton et al. [6] for variable coordination numbers, the modi- 
fied quasichemical model correctly reduces to the ideal solu- 
tion model, when the Gibbs free energy change on forming 
(A — B) -pairs AgA B equals to zero. 

In the framework of MQM(1986) (the modified quasichemi- 
cal model after Pelton and Blander [5]), the coordination num- 
bers are assumed to be constant and are selected to set the com- 
position of maximum ordering x* B and to satisfy the conditions 
s con f = at x* B when AgA B — > — °°. This results in (see also Eqs. 
(30) and (31) in Ref. [5]) 



Z B = -(x* B ln(x* B ) + (1 
Z A =Z B x* B /(l -x* B ). 



•4)ln(l-4))/(4ln(2)) 



(7) 



Similar to the examples described above, consider one mole 
of binary solution A B and assume that;tg = 1/3. In this case, 
Z B = 2.75489 and Z A = 1.37744. In the limit of immiscible 
components, AgA B — > +°°. The molar configurational entropy of 
mixing given by the curves given by MQM(1986) in this case is 
presented in Fig. 8. As shown in this figure, MQM(1986) gives 
paradoxical, negative values of the configurational entropy of 
mixing for < x B < 1/3. 

Now consider the limit of highly ordered liquid, that 
is, AgA B — > — °°. Fig. 9 shows the entropy curve given by 
MQM(1986) in this case. As seen from Fig. 9, the curve is 
entirely located in the theoretically possible range. 
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Fig. 8. Molar configurational entropy of mixing given by the modified qua- 
sichemical model (1986) in the limit of immiscible components. 
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Fig. 9. Molar configurational entropy of mixing given by the modified qua- 
sichemical model (1986) in the limit of highly ordered solution. 

2.5. Modified Quasichemical Model (2000) 

In the framework of MQM(2000) (the modified quasichem- 
ical model after Pelton et al. [6]), the configurational numbers 
are assumed to be composition-dependant. The following equa- 
tions for the coordination numbers have been suggested (see 
Eqs. (19) and (20) in Ref. [6]): 



1 
1 



l 



2«AA 



1 



n AB 



2tiaa ■ 



nAB 



Z^lriAA- 



n A B 



1 



2n B B 



nAB 



(8) 



Zf B 2n BB + n AB 



Z| A 2n BB + n AB 



Consider the limit of immiscible components (AgAB — > +°°), 
assuming that the constants Z^, Zf B , Z^ B and Zg A are chosen 



as suggested by Pelton et al. [6], namely: 
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Fig. 10 demonstrates that MQM(2000) gives paradoxical val- 
ues of the configurational entropy of mixing for the entire com- 
positional range in the limit of of immiscible components. 




Fig. 10. Molar configurational entropy of mixing given by the modified 
quasichemical model (2000) in the limit of immiscible components. 

Now consider the limit of highly ordered solution, that is, 
AgAB The entropy curve given by MQM(2000) in 

this case is presented in Fig. 11. As seen form this figure, 
MQM(2000) gives negative values of the configurational en- 
tropy of mixing in the vicinity of the composition of maximum 
ordering. 
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Fig. 11. Molar configurational entropy of mixing given by the modified 
quasichemical model (2000) in the limit of highly ordered solution. 
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2.6. Classical Quasichemical Model 

It is important to note that the classical quasichemical model 
of Fowler and Guggenheim [4] (with Z A = Z B = 1/2) correctly 
reduces to the ideal solution model when AgA B — 0. In the limit 
of immiscible components, the configurational entropy of mix- 
ing given by this model correctly reduces to the theoretically 
expected curve (s con { = 0). In all other cases, the classical qua- 
sichemical model gives values of the configurational entropy of 
mixing which are within the theoretically possible range. The 
applicability of the classical quasichemical model, however, is 
limited to those chemical solutions, in which maximum order- 
ing occurs at the composition x B = 1/2. 

2.7. Summary 

The classical associate model results in paradoxically large 
values of the configurational entropy in the limit of ideal so- 
lution, while giving theoretically possible values of the config- 
urational entropy of mixing in other two limits. The associate 
species model can give paradoxically large values in the limit 
of ideal solution and in the limit of highly ordered solution. 

The modified quasichemical model after Pelton and Blan- 
der [5] results in paradoxically small, negative values of the 
configurational entropy in the limit of immiscible components. 
The modified quasichemical model after Pelton et al. [6] gives 
paradoxically small, negative values in the limit of immiscible 
components and in the limit of highly ordered solution. 

The modified associate formalism [1] and the classical qua- 
sichemical model [4] give theoretically possible values of the 
configurational entropy in all considered limits. The classical 
quasichemical model also correctly reduces to the theoretically 
expected curve in the limit of immiscible components. The ap- 
plicability of the classical quasichemical model, however, is 
limited. 

3. Dilute Solutions 

Pelton et al. [6] considered binary solution, in which maxi- 
mum ordering occurs at the composition x B = 1/3 in the limit 
of highly ordered solution (Gibbs free energy change on form- 
ing (A — B)-pairs approaches — °°). It was demonstrated that, in 
the framework of the modified quasichemical model, chemical 
activity as of the component B in the varies in the vicinity of 
xa = as (1 — xa) even in this limiting case. 

Similar solution has also been considered in Ref. [1] in the 
framework of the modified associate formalism. As demon- 
strated, chemical activity a B varies as (1 — xa) only for finite 
values of Ag2,i and Ag\2- Here, Ag2.\ and Ag\2 are the mo- 
lar Gibbs free energy changes on forming ^B-associates and 
AZ?2-associates, respectively. In the limit of highly ordered so- 
lution (Ag2 : i — > — °°)> a B varies as (1 — xa/2). Note also that 
the interval in the vicinity of xa = 0, in which ag can be ap- 
proximated as (1 — xa), decreases with the decrease in Ag2,i- 
This should be taken into consideration, when the modified as- 
sociate formalism is used for modelling dilute solutions. 



Now consider the solution discussed in in Ref. [1] in the 
framework of the associate species model. That is, we consider 
the associate species A q , B q , A 2q /3B q /j, and A q pB 2q /3 of the 
size q. In this case, chemical activity ag is given by 



a B 



Mi 



(10) 



where x Bq is the molar fraction of the species B q . Eq. (28) in 
Ref. [1] reads 

-1 



das _ 1 i/q-i f dx A 



= —x 



dx A q 



dx Ba 



(11) 



Using the technique described in Section 6 of the previous 
paper [1], one verifies that 
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(12) 
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(f 2 + ±e 2 + feif) 

Here, t = (x Aq /x Bq ) ' /3 , £i = exp (-Ag A2q/3 B q/3 /RT) and e 2 = 

exp(-Ag As/3B2?/3 //?r), while A SA 2q/3 B q/3 and Ag AqpB2q/3 are 
the molar Gibbs free energies of formation of A 2^/3-6^/3 -species 
and A q pB 2q ^-species, respecively. 

Substituting Eq. (12) into Eq. (11) and taking the limit when 
xa — > 0, one verifies that 

3 

lim (da B /dxA) 



q 



(13) 



This demonstrates that, for any finite values of AgA 2 , 3 s , 3 and 
A?A / 3 flw3> chemical activity a B given by the associate species 
model varies as (1 — xa) only when q = 3. In the limit of highly 
ordered solution (AgA 2q/3Bq ^ — > — °°), a B varies as (1 — j^xa) in 
the vicinity of xa = 0. 

It is also important to note that, in the framework of the 
associate species model, one can consider the species A q , B q 
and A 2q pB q p only, since the set of species to be considered is 
determined by a modeller. Formally setting ^gA , 3 B 2q /3 = 
Eq. (12) reduces to 

dxA 



dx Ba 



, 2 eif (l + ieif 2 ) 
= -l-- £l f 2 + - \ . 3 V . (14) 



(f 2 + |e!f) 



Substitution of Eq. (14) into Eq. (11) and taking the limit when 

x A — > gives 

3 

lim (da B /dxA) = — — . (15) 
x A ^o 2q 

In this case, a B varies in the vicinity of xa = as (1 — j^xa) 
for any value of Ag A2q/3Bij/3 . 

The examples presented above indicate that the set of as- 
sociate species and the size of species should be carefully se- 
lected, when the associate species model is applied for dilute 
solutions. 



4. Setting the Composition of Maximum Ordering 

Modelling of multicomponent systems requires a provision 
for setting the composition of maximum ordering so as to 
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comply with the experimental data for each sub-system (bi- 
nary, ternary, quaternary and so on) individually. Associate-type 
models allow to do this in a straightforward way by including 
in each sub-system the associates, which composition coincides 
with that of maximum ordering in a particular sub-system. In 
the case of the classical associate model or the associate species 
model, a modeller can arbitrarily include the required associates 
or the associate species. In the case of the modified associate 
formalism, in which the compositions of associates are defined 
by their size, one should select the adequately large size of as- 
sociates to include the required compositions. 

In contrast to associate-type models, there are no provisions 
for setting the composition of maximum ordering for each 
ternary, quaternary and so on sub-systems individually nei- 
ther in the framework of MQM(1986) nor in the framework of 
MQM(2000). Furthermore, the composition of maximum or- 
dering can not be set independently even for binary sub-systems 
in the framework of MQM(1986). To overcome this disadvan- 
tage for binary sub-systems, Pelton et al. [6] suggested to use 
variable coordination numbers given by Eqs. (8). This modi- 
fication, however, leads to mathematical and thermodynamical 
inconsistencies of MQM(2000). These inconsistencies are de- 
scribed immediately below. 

Consider the case, in which the coordination numbers are 
given by Eqs. (8) and the constants Z^ A , Z BB , Z\ B and Z BA are 
given by Eqs. (9). Firstly, as demonstrated in Ref. [8], the total 
mole number of pairs n tp depends on the extent of the qua- 
sichemical reaction (Eq. (1) in Ref. [6]). More precisely, ntp 
decreases by 1/2 mole per each mole of (A — B) pairs formed. 
This result contradicts to the equation of the quasichemical re- 
action (see Eq. (1) in Ref. [6]), which underlies the quasichemi- 
cal model. According to Eq. (1) in Ref. [6], n t p does not change 
in the quasichemical reaction. Secondly, when the coordination 
numbers Za and Zg vary with hab, the expression for the Gibbs 
free energy of the solution given by Eq. (9) in Ref. [6]) is in 
error. Using Eqs. (1-3,19,20) in Ref. [6], one verifies that the 
correct expression in this case is 

G = ^n A g°A + ^-n B g% TAS^ + ^Ag AB . (16) 

L AA BB L 

Here, all the notations have the same meaning as in Ref. [6]. 



5. Conclusions 

The modified associate formalism has been compared with 
similar models. As demonstrated, the formalism gives theoret- 
ically possible values of the configurational entropy in all con- 
sidered limits. Being an associate-type model, the formalism 
allows to set the composition of maximum ordering for each 
sub-system of a multicomponent system individually. The cor- 
rected Gibbs free energy equation of the quasichemical model 
is also presented. 
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